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Measurements of the magnetic field penetration depth A in superconductor MoaSbT (Tc ~ 2.1 K) 
were carried out by means of muon-spin-rotation. The absolute values of A, the Ginzburg-Landau 
parameter the first Hd and the second Hc2 critical fields at T = are A(0) = 720(100) nm, 
= 55(9), ^o-H'ci(O) = 1.8(3) mT, and /io-ffc2(0) = 1.9(2) T. The zero temperature value of 
the superconducting energy gap A(0) was found to be 0.35(1) meV corresponding to the ratio 
2A(0)/fcB7c = 3.83(10). At low temperatures \~^{T) saturates and becomes constant below T ~ 
0.3Tc, in agreement with what is expected for s— wave BCS superconductors. Our results suggest 
that MoaSby is a BCS superconductor with the isotropic energy gap. 

PACS numbers: 74.70.Ad, 74.25.0p, 74.25.Ha, 76.75.+i 



Recently, the attention was devoted to MosSby. This 
compound was originally discovered more than forty 
years agoii^ and only recently was found to become a 
type-II superconductor with the transition temperature 
Tc — 2.1 K4 The properties of MosSby in a supercon- 
ducting state are rather unusual. Specific heat, resistiv- 
ity and magnetic susceptibility experiments of Candolfi et 
alA suggest that MosSby can be classified as a coexistent 
superconductor - spin-fluctuating system. As discussed 
in Rcf. factoring in the effect of spin fluctuations leads 
to renormalized values of the electron-phonon coupling 
constant and the Coulomb pseudopotential, which, be- 
ing substituted to the McMillan expression, lead to Tc 
between 1.4 K and 2.0 K. This is substantially closer to 
the experimentally observed Tc ~ 2.1 K than Tc w 10 K, 
which would be obtained without taking into account the 
effect of spin fluctuations 4 

There is currently no agreement on the symmetry of 
the order parameter of MoaSby. The recent specific heat 
experiments of Candolfi et al^ suggest that the order 
parameter in MosSby is of conventional s— wave symme- 
try. Tran et al^, based again on the results of specific 
heat measurements, have reported the presence of two 
isotropic s— wave like gaps with 2/S.i{Q) /ksTc — 4.0 and 
2A2(0)/A:bTc = 2.5 [A(0) is the zero-temperature value 
of the superconducting energy gap]. In contrast, Andreev 
reflection measurements of Dmitriev et al^^^ reveal that 
the superconducting gap is highly anisotropic. The max- 
imum to the minimum gap ratio was estimated to be 
A,„a^/A,„i„ ~ 40 and s + wave symmetry of the or- 
der parameter was proposed in a qualitative analysis^i^ 

The symmetry of the superconducting order parameter 
can be probed by measurements of the magnetic pene- 
tration depth A. A fully gaped, isotropic pairing state 
produces a thermally activated behavior leading to an 
almost constant value of the superfluid density ps oc A^^ 
for T < O.STc'^'^ Presence of line nodes in the gap leads 
to a continuum of low laying excitations, which result in 
a linear A"^(T) at low tcmperatureSfi^ii^ In two-gap su- 



perconductors with highly different gap to Tc ratios the 
inflection point in A~^(T) is generally present F^^ii^ii^ 

In this paper, we report the study of the magnetic field 
penetration depth in superconductor MoaSby by means 
of muon-spin rotation. Measurements were performed 
down to 20 mK in a series of fields ranging from 0.02 T 
to 0.2 T. Our results are well explained assuming conven- 
tional superconductivity with the isotropic energy gap in 
agreement with the recent specific heat experiments of 
Candolfi et al^ The zero-temperature value of the gap 
was found to be A(0) = 0.35(1) mcV corresponding to 
the ratio 2A(0)/A:bTc = 3.83(10). 

The MosSbr single-crystal samples were grown 
through peritectic reaction between Mo metal and liquid 
Sb.— The transverse field muon-spin rotation (TF-/iSR) 
experiments were performed at the 7rM3 beam line at 
Paul Scherrer Institute (Villigen, Switzerland). For our 
experiments the ensemble of some sub- and millimeter 
size single crystals were mounted onto the silver plate 
to cover the area of approximately 50 mm^. The silver 
sample holder was used because it gives a nonrclaxing 
muon signal and, hence, only contributes as temperature 
independent constant background. The crystals were ori- 
ented so that the magnetic field was preferably applied 
along the 001 crystallographic direction. The MosSby 
samples were field cooled from above Tc down to ~ 20 mK 
in magnetic fields ranging from 20 mT to 0.2 T. A full 
temperature scan (from 20 mK up to 2.5 K) was per- 
formed in a field of p^H = 0.02 T. 

Figure [T] shows the muon-spin precession signals in 
poH = 0.02 T above (T = 2.2 K) and below (T = 0.05 K) 
the superconducting transition temperature of MoaSby. 
The difference in the relaxation rate above and below 
Tc is due to the well-known fact that type-II supercon- 
ductors exhibit a fiux-line lattice leading to spatial in- 
homogeneity of the magnetic induction. As is shown by 
Brandt the second moment of this inhomogcneous 
field distribution is related to the magnetic field penetra- 
tion depth A in terms of (AB^) cx a^c ^ 
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FIG. 1: (Color online) Transverse-field muon-spin precession 
signals from MoaSbr obtained in hqH = 0.02 T above the 
transition temperature {T = 2.2 K - black squares) and after 
field-cooling the sample below Tc {T — 0.05 K - red circles). 
The solid lines correspond to the fit by means of Eq. [1] For 
visualization purpose the apparent precession frequencies are 
modified from the actual precession frequencies by the use of 
a rotating reference frame. 



The analysis of TF-/iSR data was carried out in the 
time-domain by using the following functional formJ^ 



A{t) ^ yloexp 



cos (7^ i3„t t -f ( 



-Af,g C0s(7^Bej;t t 



(1) 



where the first term denotes the contribution from the 
MosSbr sample and the second term is the background 
contribution from the Ag sample holder. ^0 and A\,g are 
the initial asymmetries arising from the sample and the 
background, -Bmt and -Bext are the internal filed inside 
the sample and the applied external field seen in the Ag 
backing plate, 7^ = 27r x 135.5342 MHz/T is the muon 
gyromagnetic ratio, and is the initial phase. Osc and 
Onra are the muon-spin relaxation rates caused by the 
nuclear moments and the additional component appear- 
ing below Tc due to nonuniform field distribution in the 
superconductor in the mixed state. 

The analysis was carried out as follows. First, the 
muon-time spectra were fitted by means of Eq. (HI) with 
all the parameters free. Then, the ratio Ao/Atg = 0.975 
was obtained as the mean value in the temperature in- 
terval of 0.02 K<T < 1.5 K. After that, by keeping the 
ratio Abg/Ao fixed, the nuclear moment contribution anm 
was estimated by analyzing the experimental data above 
2.1 K, (i.e. at T > Tc, where cr^c = 0). The final analysis 
was made with Ao/Abg, anm and B^xt as fixed, and Bint 
and (T,5c as free parameters. This allows to reduce the 
number of the fitting parameters and, as a consequence, 
to increase the accuracy in determination of ctsc- 

Within BCS theory the zero-temperature magnetic 
penetration depth [A(0)] of isotropic superconductor does 
not depend on the magnetic field. Contrary, the nonlocal 
and the nonlinear response of the superconductor with 



nodes in the gap to the magnetic field, as well as, the 
faster suppression of the contribution of the smaller gap 
to the total superfiuid density in a case of two-gap super- 
conductor leads to the fact that A(0), evaluated from /iSR 
experiments, is magnetic field dependent and it increases 
with increasing fieldj ^°i^^ Such behavior was observed in 
various hole-doped cupratesr^i22ii^ and the double-gap 



NbSe- 



22.24 



0.20 



and MgB2 superconductors 



25.26 



0.15 




0.10 



100 150 
6(mT) 

FIG. 2: (Color online) Field dependence of the fiSH depolar- 
ization rate CTsc measured at T = 20 mK in MosSby sample. 
The solid line corresponds to a fit of Eq. ([2]) to the experimen- 
tal data with the parameters shown in the figure. The inset 
shows the dependence of normalized to its value at B = 
on the reduced magnetic field b — B/Bc2 for MosSbr, the two- 
gap superconductor MgB2 from Ref. [2^, and the hole-doped 
cuprate superconductor YBa2Cu3 06.95 from Ref. [2^. 

Consideration of the ideal vortex lattice (VL) of 
an isotropic s— wave superconductor within Ginsburg- 
Landau approach leads to the following expression for 
the magnetic field dependence of the second moment of 
the magnetic field distribution>i^ 



<7sc[l-^s ^] = a X (1 - B/Bc2) 

1 + 1.21(1- y^B/Bc2^^ 



A ^ [nm] 



(2) 



Here B is the magnetic induction, which for applied fields 
in the region Hd < Happ < iJc2 is i? ~ fJ-oHapp, a 
is the coefficient depending on the symmetry of the VL 
(a = 4.83 X lO'' for triangular VlJ^d^ and, as is shown 
below, a = 5.07 x 10^ for rectangular VL), Hd is the 
first critical field, and Bc2 = fJ-oHc2 is the upper criti- 
cal field. Equation ([2]) accounts for reduction of age due 
to stronger overlapping of vortices by their cores with 
increasing field. According to calculations of Brandtfi^ 
Eq. ^ describes with less than 5% error the field vari- 
ation of (Tsc for an ideal triangular vortex lattice and it 
holds for type-II superconductors with the value of the 
Ginzburg-Landau parameter K = A/^>5(^is the coher- 
ence length) in the range of fields 0.25/k^ '^ < B/Bc2 < 1- 
Satisfactory fit of Eq. ([2]) to the experimental data 
would suggest that there is no significant change of the 
penetration depth in the range of magnetic fields of 
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the experiment. We performed an analysis of the mag- 
netic field dependence of the muon depolarization rate 
Use measured at T = 20 mK (see Fig. which is a 
good approximation of OsJyF = 0, -B). The solid line in 
Fig. [2] corresponds to the fit of experimental (Jsc{B) by 
means of Eq. Q. The fit yields A(0) = 720(25) nm and 
Bc2{0) — 1-9(2) T. A good agreement between the ex- 
periment and the theory (see Fig. ^ allows to conclude 
that asc{T = 0,5) can be consistently fitted within the 
behavior expected for conventional isotropic supercon- 
ductor. For comparison in the inset of Fig. [2] we plot 
the dependence of the inverse squared penetration depth 
normalized to its value at B = on the reduced field 
b = B/Bc2 for MoaSby studied here, the two-gap su- 
perconductor MgB2 from Ref. [2^, and the hole-doped 
cuprate superconductor YBa2Cu3 06.95 from Ref. |22. It 
is seen that in a case of MgB2 and YBa2Cu3 06.95 A~^, 
evaluated from muon experiments, is field dependent and 
decreases rather strongly with increasing field. 
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FIG. 3: (Color online) (a) Theoretical P{B) profiles (A = 
720 nm, C = 13 nm, and ^loH = 0.02 T) calculated for 
(jg = 0.0006, 0.05, 0.1, and 0.15 /is"\ (b) Dependence of 
the product dsc x A~^ on the reduced field 6 = B /Bc2 for the 
triangular (circles) and the rectangular (stars) VL's. Calcula- 
tions were made by using the numerical approach of Brandlii^ 
for K = 100. The solid lines are obtained by means of Eq. ((2]) 
with a = 4.83 • 10* for triangular and a = 5.07 • 10* for rect- 
angular VL's, respectively. The insets are the contour plots 
of field variation within triangular and rectangular VL's. 

The absolute error in A(0), estimated from the fit of 
asc{B) by means of Eq is ~ 3.5%. There are, how- 
ever, other sources of uncertainties which could also af- 
fect A(0) and, consequently, increase the error: (i) The 
use of Eq. ^ to extract A^^ from asc{B) requires that 
the VL is well ordered and it is triangular, as opposed 
to disordered or rectangular VL. (ii) The fitting function 
Eq. (HI) assumes that the contribution from the sample 
is described by a single line of Gaussian shape as op- 
posed to asymmetric magnetic field distribution P{B) 
generally observed in good quality single crystals 32, (iii) 
The background signal due to the Ag backing plate was 
assumed to have a temperature-independent frequency 
and is non-relaxing at all T. However, when the sample 
goes superconducting it excludes some field, even in the 
mixed state. These small fields are very inhomogeneous 



and could cause as a shift of the background signal to 
higher fields and an increase of a slow relaxation. In the 
following we are going to discuss in detail the sources (i) 
and (ii). We believe that the source (iii) does not play 
a substantial role here since even at a T = 20 mK Bint 
is only 0.16 mT smaller than the applied field and Bg^t 
is constant in the whole temperature range (from 20 mK 
up to Tc). It should be mentioned, however, that experi- 
ence with more robust superconductors has shown that 
even though the background field does not usually change 
perceptibly, the excluded flux can affect the background 
relaxation rate. 

To account for possible random deviations of the flux 
core positions from their ideal ones (VL disorder) and for 
broadening of /iSR spectra due to nuclear depolarization, 
the field distribution of an ideal VL [Pid{B)] is, gener- 
ally, convoluted with a Gaussian distribution in terms 

Qf .27.28.29 



P{B) = 
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and ub is the contribution to 



the Gaussian broadening of Pid{B) due to VL disorder. 
The theoretical P{B) profiles for various Ug's are shown 
in Fig. [3] (a). It is obvious that both, the nuclear mo- 
ment contribution and the VL disorder broaden the P{B) 
profiles, thus requiring that the total second moment of 
/iSR line needs to be the sum of three components)^ 
+ <^nm + "'b- This implies that neglecting the VL 
disorder leads to overestimation of cTgc and, as a con- 
sequence, to underestimation of A. On the other hand, 
asc, obtained from the fit of asymmetric P{B) line by us- 
ing symmetric Gaussian function [see Eq. ([T])], becomes 
underestimated.— In a case of MosSby the main con- 
tribution to Ug comes, most probably, from the nuclear 
moment term cr„m ~ 0.178 /is~^, which is comparable 
with Use in the whole field range, rather than from ctb, 
which for good quality single crystals corresponds, typi- 
cally, to 10 — 20%, of cTjc (see, e.g., Ref. [2^ and references 
therein). Such big anm also leads to the fact that the 
shape of P{B) profile is close to the Gaussian one as is 
demonstrated by the solid blue line in Fig. [3] 

Dependences of <Jsc on the reduced field b — B /Bc2 for 
triangular and rectangular VL's, obtained by using nu- 
merical calculations of Brandt are shown in Figl3] (b). 
Solid hues correspond to Eq. ^ with a = 4.83 x 10* (red 
line) and 5.07 x 10* (blue line). As is seen, in a case of 
rectangular VL <Jsc{b) can still be satisfactory described 
by means of Eq. ([2|) with the only 5% bigger value of the 
coefficient a. This implies that the uncertainty related 
to different VL symmetries would result in ~ 2.5% addi- 
tional error in A. To summarize, by taking into account 
the above presented arguments we believe that the true 
overall uncertainty in the absolute value of A(0) is about 
10-15 %, and A(0) = 720(100) nm. 

The zero-temperature value of the superconducting 
coherence length ^(0) may be estimated from Bc2{0) 
as ^(0) = [$o/27^Bc2(0)]"■^ which results in ^(0) = 
13(1) nm (<f>o is the magnetic fiux quantum). Using the 
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values of A(0) and ^(0), the zero-temperature value of 
the Ginzburg-Landau parameter is k(0) = A(0)/^(0) = 
55(10). The value of the first critical field can also be cal- 
culated by means of Eq. (4) from Ref. [H as /ioi?ci(0) = 
1.8(3) mT. Note that the zero-temperature values of A(0), 
^(0); k{0), and the first and the second critical fields are 
in agreement with those reported in the literaturci^i^i^'^ 
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FIG. 4: (Color online) Temperature dependence of 
of MoaSby reconstructed from asc{T) measured in (IqH — 
0.02 T. The black solid and the dash-dotted blue lines are the 
fits by means of clean and dirty weak-coupling BCS models. 
The blue dashed line represents A~^(r) calculated by using 
two-gap alpha mode l^^'^^ with the parameters from Ref. [g. 
See text for details. 



As a next step we are going to analyze the temper- 
ature dependence of the magnetic penetration depth. 
A~^(T), presented in Fig. SI was reconstructed from 
(Tsc{T) measured in a field of 0.02 T by using Eq. 
BcaiT), needed for such reconstruction, was assumed to 
follow the Werthamer-Helfand-Hohenberg prediction^ in 
agreement with the results of Ref. H. It is worth to 
mention here that the temperature dependence of A~^ 
is much less affected by the above described uncertain- 
tics than the absolute A(0) value. Indeed, both, as^^^ 
and the relaxation related to the excluded flux increase 
with decreasing T in much the same way as CTsc, which 
causes mainly a correction to the coefficient a in Eq. ([2]), 
while (T„m is temperature independent. 

Figure [3] implies that below 0.7 K A^^ is tempera- 
ture independent as expected for the superconductor with 
fully gaped states. The experimental X~^{T) depen- 
dence was analyzed within the dirty- and the clean-limit 
approaches. In the dirty-limit theory A~^(T) has the 
formi^ 
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Here / = [1 -|-exp(i?/fcBT)] ^ is the Fermi function. The 
temperature dependence of the gap was approximated 
by A(r) = A(0)tanh{1.82[1.018(r,/T - As 
is seen, both, the dirty- and the clean-limit approaches, 
describe the experimental data reasonably well. The fit- 
ted curves are almost undistinguishable from the each 
other. The results of the fits are =2.11(2) K, 
A(0) =0.35(1) meV and 2.12(2) K, 0.39(1) meV for 
the dirty- and the clean-limit BCS model, respectively. 
The corresponding gap to Tc ratios were found to be 
[2^{Q)/kBT,]d„ty = 3.83(10) and [2A(0)/A:bTJ 

clean — 

4.27(11). There are few reasons why, we believe, MosSby 
studied here is in the dirty limit: (i) The maximum 
gap value obtained in Andreev reflection experiments 
was found to be A(0) ~ 0.32 meV which is close to 
A(0) =0.35(1) meV observed within our dirty limit cal- 
culations, (ii) The residual resistivity ratio measured 
on similar single crystals was found to be rather low 
P300k/pok = 1.44 Note that the isostructural com- 
pound RusSnr has Psook/pok = 144;^ which is bigger 
by more than 2 orders of magnitude, (iii) A reason- 
able assumption about Fermi velocity vp ~ 10^ m/s can 
be made within the free electron approximation by tak- 
ing Ep — 6.5 eV<^ This allows us to readily estimate 
the BCS-Pippard coherence length = hv p / it ~ 
600 nm which is approximately 50 times bigger than 
^ = 13(1) nm obtained experimentally. Considering that 
=ec7^+ri {I is the mean free path) we obtain that 
the coherence length ^ in MoaSby is limited by I and, 
correspondingly, Z ~ ^ ~ 13 nm. 

Now we are going to comment shortly the results 
of the speciflc heat experiments of Tran et al^ It was 
shown that the analysis of the electronic speciflc heat 
data within the framework of the phenomenological 
alpha model suggests the presence of two gaps with 
2Ai(0)/A:bTc = 4.0 and 2A2(0)/fcBrc = 2.5 and the rel- 
ative weight of the bigger gap of 0.7fS, The \~'^{T) curve 
simulated by using these parameters and A(0) = 716 nm, 
is shown in Fig.|4]by blue dashed line. Both contributions 
were assumed to be in the dirty limit [see Eq. @] and the 
similar alpha model, but adapted for calculation of the 
superfluid density, was usedj ^'^d^ The agreement between 
the simulated curve and the experimental data is rather 
good. We should emphasize, however, that in two-gap su- 
perconductor the contribution of the smaller gap to the 
total superfluid density decreases very fast with increas- 
ing held; thus leading to strong suppression of is, 25^26 
This is inconsistent with the data presented in the inset 
of Fig. [5] revealing that in MoaSby X~'^{B) = const. 

To conclude, the superconductor MosSbj (Tc ~ 2.1 K) 
was studied by means of muon-spin rotation. The main 
results are: (i) The absolute values of the magnetic fleld 
penetration depth A, the Ginzburg-Landau parameter k, 
and the first Hd and the second Hc2 critical fields at 
T = were found to be A(0) = 720(100) nm, k(0) = 
55(8), noHciiO) = 1.8(3) mT, and noHc2iO) = 1-9(2) T. 
(ii) Over the whole temperature range (from Tc down to 
20 mK) the temperature dependence of A~^ is consistent 
with what is expected for a single-gap s— wave BCS su- 
perconductor, (in) The ratio 2A{0) /ksTc = 3.83(10) was 
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found, suggesting that MosSby superconductor is in the 
intermediate-couphng regime, (iv) The magnetic pene- 
tration depth A is field independent, in agreement with 
what is expected for a superconductor with an isotropic 
energy gap. (v) The value of the electronic mean-free 
path was estimated to be Z ~ 13 nm. This relatively 
short value suggests that strong scattering processes play 
an important role in the electronic properties of MoaSby. 
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